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Abstract 



It has been proved that the spanning tree from a given network has the optimal synchronizability, which means the index R - A^l 
reaches the minimum 1 . Although the optimal synchronizability is corresponding to the minimal critical overall coupling strength 
Qs^ to reach synchronization, it does not guarantee a shorter converging time from disorder initial configuration to synchronized state. 
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In this letter, we find that it is the depth of the tree that affects the converging time. In addition, we present a simple and universal 
|way to get such an effective oriented tree in a given network to reduce the converging time significantly by minimizing the depth of 
the tree. The shortest spanning tree has both the maximal synchronizability and efficiency. 

'Key words: Synchronization, Spanning tree, Kuramoto model 
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1. Introduction 

' The synchronization is an universal phenomenon emerged 
by a population of dynamically interacting units. It plays 
Sn important role from physics to biology and has attracted 
much attention for hundreds of years. Thus, there are a great 
deal of relative researches based on this topic. With the un- 
derstanding of relations between network topology and the 
synchronizability[l-5], scientists have proposed many meth- 
ods to enhance synchronizability of the network[6-15]. Some 
t)f them tried to modify the topology of the network to en- 
hance the synchronization[6-10] while others by just modifying 
the coupling weight of each edge while keeping the topology 
unchanged[ 11 - 1 6] . 

In these papers, the synchronization is always measured by 
the eigenvalues of Laplacian matrix as /? = I A^iill -19}. The 
smaller the R is, the better the synchronizability will be. In re- 
cent years, a lot of works focus on how to enhance the synchro- 
nizability by distributing the weight to the edges according to 
the structural characteristics of the nodes and edges. For exam- 
ple, distributing the weight by the degree and betweeness can 
pharply enhance the network synchronizability, these methods 
can reduce index R io a. small value[l 1-14]. For the growing 
scale-free network, some researchers took the age into consid- 
eration and reduced the index R to an even smaller value close 
to the minimum 1[15-16]. In the Ref[20], Nishikawa and Mot- 
ter gave the weight distributing an extreme way by imposing 
the weight of some edges to 0. The process can be regarded 
as cutting off the edges which are disadvantage to the synchro- 
nization. Finally, they can get an oriented tree with normalized 
input strength and no directed loops. Moreover, they proved 
that the R of the tree is 1, meaning that the tree has maximal 
synchronizability. The index R is corresponding to the critical 



overall coupling strength. When R reaches its minimum 1, the 
synchronized states are stable for the widest possible range of 
the parameter representing the overall coupling strength. 

When investigating the synchronization of a network, we al- 
ways consider both of the critical overall coupling strength to 
synchronized the whole network and the converging time. The 
converging time can be regarded as the efficiency of a network. 
However, better synchronizability does not guarantee a shorter 
converging time from disorder initial configuration to synchro- 
nized state. Actually, in the Ref[20], Nishikawa and Motter 
found that the synchronizing process may take longer time in 
the optimal network with R - \. li leads us to an interest- 
ing problem that what is the factor that affects the converging 
time. In this paper, we will give a clear answer to this problem. 
In addition, a simple and universal method is presented to find 
the optimal spanning tree with maximal synchronizability and 
efficiency. This kind of tree is the optimal structure for synchro- 
nization that the scientists try to find from any given network. 



2. Result 

2.1. The factor affecting the converging time 

In a dynamical network, each node represents an oscillator 
and the edges represent the couplings between the nodes. For a 
network of linearly coupled identical oscillators, the dynam- 
ical equation of each oscillator can be written as 
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It has been proved that the synchronizability of an oriented 
tree of the network has reaches its maximum, that is the index 
R = 1 and the synchronized states are stable for the widest pos- 
sible range of overall coupling strength. Generally, this kind 
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Figure 1: The converging time to synchronization under trees with different 
depth. All the trees are of the same size = 10. The only different between 
them is the depth of the tree. The □ stands for the overall coupling strength tr is 
5, the O stands for cr = 10 and the A stands for cr = 20 .The result is averaged 
by 50 times. 

of oriented tree has a root with no input. It works as the mas- 
ter oscillator and affects the oscillators in the hierarchical level 
below without any feedback. Then, the next lower level oscil- 
lators will get synchronized and so on until the whole network 
reaches complete synchronization. Hence, the hierarchy num- 
ber is very important for it determines the converging time[21]. 
Clearly, in this kind of spanning tree, the synchronized process 
in any branch is independent with each other. So the oscillator 
in the lowest hierarchical level would reach the synchronized 
state at last. From this point of view, we figure that when given 
the oscillator model and the overall coupling strength, the con- 
verging time is only determined by the depth of the tree. 

In order to validate the assumption, we put the Kuramoto 
model to each node of the network to make numerical simula- 
tion. The Kuramoto model is a classical model to investigate 
the phase synchronization phenomenon[22-24]. The coupled 
Kuramoto model in the network can be written as 

Oi = oji + crGijSiniej - 0i), i = 1,2 ,N (2) 

The collective dynamics of the whole population is measured 
by the macroscopic complex order parameter, 

r(0e'*« = ^2e''^« (3) 

7=1 

Where the r(t) ^ 1 and r(t) ^ describe the limits in which 
all oscillators are either phase locked or move incoherently, re- 
spectively. In this paper, we define the time T as the converging 
time, where for any t > T, r{t) - r{t - dt) < s all along and the 
s = 10-^ 

To investigate which affects the converging time of the net- 
work, we adopt the trees with 10 nodes. When the depth is 1, 
the tree has only one root, all the other nodes which only receive 
input from the root are located in the second hierarchical level. 



Figure 2: The converging time to synchronization under trees with different 
size. The overall coupling strength <t is kept the same as 10. The A stands for 
the depth of the tree equals 1 , the O stands for depth = 5 and the □ stands for 
depth = 9 .The result is averaged by 50 times. 

WTien the depth of the tree is 9, the tree is just a chain connect- 
ing all the nodes. When the depth of the tree is from 2 to 8, the 
structures are more complicated. We compare the converging 
time of the trees with different depths. Since the trees with the 
same depth may have different structures, we use the average T 
to represent the converging time. From Fig.l, it is clear that the 
deeper the tree is, the more converging time is needed to reach 
synchronized state. The effect of the depth to converging time 
is Unear roughly. 

Furthermore, given the depth, the converging time will not be 
affected by the size of the tree. It is reasonable because every 
branch of the tree is independent, and the longest branch will 
reach the synchronized state at last no matter how many the 
branches are. It can be seen from the compare between the 
trees from 10 nodes to 100 nodes from Fig. 2. In the simulation, 
we keep the two kind of trees with the same depth such as 1, 5 
and 9. 

So it is the depth of the tree that actually affects the con- 
verging time of the synchronization process. Additional, when 
given the depth of a tree, with specific oscillator model and 
overall coupUng strength cr, the converging time is almost a 
constant. Since the initial state of each oscillator is given ran- 
domly, there may be some fluctuation in the constant, which 
can be described by the error. Thus, if we want to find a best 
way to distribute the weight to enhance the synchronization, 
we just have to look for the spanning tree with minimal depth. 
This kind of tree has minimal R value and shortest converging 
time, which means that it has maximal synchronizability and 
efficiency. 

2.2. Center of the network and the shortest spanning tree 

To get such a tree is to create a shortest oriented tree from 
a undirected network. Each tree has a root, the root adopting 
is very important because it determines the depth of the tree. 
However, the root is not simply the node with the largest degree. 
For in some case, the spanning tree created by this root may be 
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Figure 3: Given a network with 10 nodes, the left tree is created by adopting 
the node with largest degree as the root, and the right one is the shortest tree 
based on this topology. Obviously, the center is not simply the node with the 
largest degree. 

not the shortest. A typical example is shown in Fig. 3. Here, we 
call the node to create the shortest spanning tree as the center 
of the network. 

We use the signaling process to find the center For a network 
with nodes, every node is assumed to be a system which can 
send, receive, and record signals. A node can only affect its 
neighbors, which will affect their neighbors in the same way. 
Finally, each node will affect the whole network. At the begin- 
ning, we set every node as the source and give each of them one 
unit of signal, the signaling process is independent. For each 
step, the node will transmit the signal to the neighbors. For the 
signaling process is independent, several signal can be trans- 
mitted in a node simultaneously. After m time steps, there must 
be a signal which covers all the nodes in the network with the 
fewest steps. So corresponding source node is the center of the 
network which is the root of the shortest spanning tree. 

Actually, the above signaling process could be described by 
a simple but clear mathematical mechanism [25]. Suppose we 
have a network with nodes, it can be represented mathemat- 
ically by an adjacency matrix A with elements A,y equals to 1 
if there is an edge from / to j and otherwise. And the / is a 
matrix where all the diagonal elements are 1 and others are 0. 
Then the column / of matrix V = (I + A)™ will represent the 
effect of source node / to the whole network in m steps. So we 
can get a A^-dimensional vector that records each node's signal 
quantity which represents the effect of the source node. If all 
the elements of a column are nonzero, the signal of the source 
node / has affected the whole network. For each step, all the 
columns will be updated. To find the center, we should simply 
find which column reaches totally nonzero at first. 

The shortest spanning tree can be created by the signaling 
process too. Suppose node / is the center of the network, we 
mark ; as the used node. The second level includes all the un- 
used nodes connected to the node /, then mark all these nodes 
and the edges from / to them as used. To create the third level, 
we consider each node in the second level as sub-center, and 
the process is the same as the center /. Specifically, the edges 
between the nodes in the same level should be left unmarked. 
Moreover, if a node k is marked by one of the node in the higher 




Figure 4: The left tree is a shortest spanning tree created from a scale-free 
network, which is from BA model with average degree 6. The right one is from 
a random network with the same average degree as the scale-free network. It 
clear that the depth of such effective directed trees is related to the topology of 
the original networks. 

level, although the node k is connected to another node j in 
the higher level, the edge between k and j should be left un- 
marked too. After all the nodes are marked, we can get the 
shortest spanning tree. The nodes in the tree are connected by 
the marked edge and the direction of the edge is from the higher 
level to the lower level. This spanning tree (i) embeds a directed 
spanning tree, (ii) has no directed loop, and (iii) has normalized 
input strengths as the Ref[20]. So the index R of the tree equals 
to 1 . Additionally, For its depth is minimal in all spanning trees 
from the given network, so the synchronization converging time 
will be the shortest. 

The topology of the original network is related to the depth 
of such effective directed trees. For instance, the center of the 
scale-free network usually has big degree. This kind of center 
can reduce the depth of the tree significantly. On the contrary, 
the homogenous random network does not have such kind of 
center, so the depth of the spanning tree from these networks 
will always be longer than that from the scale-free network as 
shown in Fig .4. 

3. Conclusion 

In the former works, to enhance the synchronization is al- 
ways to reduce the index R. That is to reduce the overall cou- 
pling strength. However, the converging time is also an impor- 
tant factor In this paper, we find that the depth of the tree is 
the only factor that effects the converging time when given the 
oscillator model and the overall coupling strength cr. Addition- 
ally, we present a simple way to obtain the shortest spanning 
tree with maximal synchronizability and efficiency from any 
given network. 

In order to enhance the synchronizability of a network, the 
coupling strength of each edge is always be scaled by some 
structural characteristics such as degree, betweeness and so 
forth. The purpose of this process is to reduce the effect of 
the edges which is disadvantage for synchronization. To make 
the process to an extreme way, creating a spanning tree from 
a network is to cut off some edges in the network, which can 
also be considered as imposing the weight of these edges to 0. 
All the spanning trees from the network have the same R - 1, 
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which means the critical overall couple strength of the tree is 
minimal. Hence, the tree creating is a coupling strength scaling 
process. 

Moreover, compared with any spanning tree from a given 
network, the shortest spanning tree can reduce the converging 
time significantly. It makes the tree have maximal efficiency. 
In fact, comparing with the former works[15-16], the way pre- 
sented by us is universal. It is valid not only in the heteroge- 
neous networks, but also in the homogenous networks. 
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